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Abstract 

We present an exact solution for a catalytically-activated annihilation A + A — » 
reaction taking place on a one-dimensional chain in which some segments (placed at 
random, with mean concentration p) possess special, catalytic properties. Annihila- 
tion reaction takes place, as soon as any two A particles land from the reservoir onto 
two vacant sites at the extremities of the catalytic segment, or when any A particle 
lands onto a vacant site on a catalytic segment while the site at the other extremity 
of this segment is already occupied by another A particle. We find that the disorder- 
average pressure p(i uen ) per site of such a chain is given by p(i uen ) = p(lan) _|_ (3~ 1 F, 
where = ln(l + z) is the Lan gmuir adsorption pressure, (z being the ac- 

tivity and - the temperature), while [3~ 1 F is the reaction-induced contribution, 
which can be expressed, under appropriate change of notations, as the Lyapunov 
exponent for the product of 2 x 2 random matrices, obtained exactly by Derrida and 
Hilhorst (J. Phys. A 16, 2641 (1983)). Explicit asymptotic formulae for the particle 
mean density and the compressibility are also presented. 

PACS numbers: 82.65.+r; 64.60.Cn; 68.43.De 



1 Introduction. 



Catalytically-activated reactions (CARs), i.e. reactions between chemically inactive molecules 
which recombine only when some third substance - the catalytic substrate - is present, are 
widespread in nature [1,2]. Recently, such reactions have attracted a considerable atten- 
tion following an early observation [3] of remarkable non-mean-field behavior exhibited 
by a specific reaction - the CO-oxidation in the presence of metal surfaces with catalytic 
properties [1,2]. An extensive analysis of this CAR has substantiated the emergence of an 
essentially different behavior as compared to the predictions of the classical, formal-kinetics 
scheme and have shown that under certain conditions such collective phenomena as phase 
transitions or the formation of bifurcation patterns may take place [3] . Prior to these works 
on catalytic systems, anomalous behavior was amply demonstrated in other schemes [4], 
involving reactions on contact between two particles at any point of the reaction volume 
(i.e., the "completely" catalytic sysems). It was realized [4] that the departure from the 
text-book, formal-kinetic predictions is due to many-particle effects, associated with fluc- 
tuations in the spatial distribution of the reacting species. This suggests that similarly to 
such " completely" catalytic reaction schemes, the behavior of the CARs may be influenced 
by many-particle effects. 

Apart from the many-particle effects, behavior of the CARs in practically involved 
systems might be affected by the very structure of the catalytic substrate, which is often 
not well-defined geometrically, but must be viewed as an assembly of mobile or localized 
catalytic sites or islands, whose spatial distribution is complex [1]. Metallic catalysts, 
for instance, are often disordered compact aggregates, the building blocks of which are 
imperfect crystallites with broken faces, kinks and steps. Another example is furnished by 
porous materials with convoluted surfaces, such as, e.g., silica, alumina or carbons. Here 
the effective catalytic substrate is also only a portion of the total surface area because of 
the selective participation of different surface sites to the reaction. Finally, for liquid-phase 
CARs the catalyst can consist of active groups attached to polymer chains in solution. 

Such complex morphologies render the theoretical analysis difficult. As yet, only em- 
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pirical approaches have been used to account for the impact of the geometrical complexity 
on the behavior of the CARs, based mostly on heuristic concepts of effective reaction order 
or on phenomenological generalizations of the formal-kinetic "law of mass action" (see, e.g. 
Refs.fl] and [2] for more details). In this regard, analytical solutions of even somewhat 
idealized or simplified models, such as, for instance, the ones proposed in Refs.[3], are 
already highly desirable since such studies may provide an understanding of the effects of 
different factors on the properties of the CARs. 
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Figure 1: One-dimensional lattice of adsorption sites in contact with a reservoir. Filled circles 
denote hard-core A particles. Thick black lines denote the segments with catalytic properties, 
(a) denotes a "forbidden" particle configuration, which corresponds to immediate reaction, (b) 
depicts the situation in which two neighboring A particles may harmlessly coexist. 



In this paper we study a catalytically-activated annihilation A + A — > reaction in 
a simple, one-dimensional model with random distribution of the catalyst, appropriate to 
the just mentioned situation with the catalytically-activated reactions on polymer chains. 
We present here an exact solution for this model with quenched random distribution of the 
catalyst and show that despite its apparent simplicity it exhibits an interesting non-trivial 
behavior. We note finally that kinetics of A + A — > reactions involving diffusive A par- 
ticles which react upon encounters on randomly placed catalytic sites has been discussed 
already in Refs. [5-7] and [8] , and a rather surprising behavior has been found, especially in 
low-dimensional systems. Additionally, steady-state properties of A + A — > reactions be- 
tween immobile A particles with long-range reaction probabilities in systems with external 
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particles input have been presented in Refs.[9] and [10] and revealed non-trivial ordering 
phenomena with anomalous input intensity dependence of the mean particle density, which 
agrees with experimental observations [11]. 



2 The model. 



Consider a one-dimensional regular lattice of unit spacing comprising iV adsorption sites. 
The lattice is in contact with a reservoir of identical, non-interacting hard-core A particles 
(see, Fig.l) - a vapor phase, which is steadily maintained at a constant pressure. 

The A particles from the vapor phase can adsorb onto vacant adsorption sites and 
desorb back to the reservoir. The occupation of the "i"-th adsorption site is described by 
the Boolean variable rii, such that 



Suppose next that some of the segments - intervals between neighboring adsorption sites 
possess "catalytic" properties (thick black lines in Fig.l) in the sense that they induce an 
immediate reaction A + A — > 0, as soon as two A particles land onto two vacant sites at 
the extremities of the catalytic segment, or an A particle lands onto a vacant site at one 
extremety of the catalytic segment while the site at the other extremity of this segment 
is already occupied by another A particle. Two reacted A particles instantaneously leave 
the lattice (desorb back to the reservoir). Any two A particle adsorbed at extremities of a 
non-catalytic segment harmlessly coexist. 

To specify the positions of the catalytic segments, we introduce a Boolean variable Q, 
so that Co — (n = and 



In what follows we suppose that Q are independent, identically distributed quenched ran- 




0, otherwise. 



1, if the "i"-th site is occupied, 




1, if the i-t interval is catalytic, i = 1,2, . . . ,N — 1, 
0, otherwise. 
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dom variables with distribution 



p( C )=p5(C-l) + (l-p)5(C). (!) 

Now, for a given distribution of the catalytic segments, the partition function Z N (() of 
the system under study can be written as follows: 

Zn(C) = E ^ II i 1 - C* » (2) 

{nj i=l 

where the summation extends over all possible configurations {n^}, while 

z = exp(/3/x) (3) 

is the activity and /i - the chemical potential. Note that Z N (() in Eq.(2) is a functional of 
the configuration ( = 

It is worth-while to remark that Z N (() can also be thought of as a one-dimensional 
version of models describing adsorption of hard-molecules [12-19], i.e. adsorption limited 
by the "kinetic" constraint that any two molecules can neither occupy the same site nor 
appear on the neighboring sites. The most celebrated examples of such models are fur- 
nished by the so-called "hard- squares" model [12-16], or by the " hard- hexagons" model 
first solved exactly by Baxter [18]. In our case of the CARs on random catalytic substrates 
the nearest-neighbor exclusion constraint is introduced only locally, at some specified, ran- 
domly distributed intervals. Such locally frustrated models of random reaction/adsorption 
thus represent a natural and meaningful generalization of the well-studied exclusion mod- 
els over systems with disorder. Of course, in this context two-dimensional situations are 
of most interest, but nonetheless it might be instructive to find examples of such models 
which can be solved exactly in one dimension. 

Our main goal here is to calculate the disorder- average pressure per site: 

p<, " e "S^( ln(ZK(c)) )<- (4) 

where the angle brackets with the subscript ( denote averaging over all possible configu- 
rations {Ci}- Once p(9" en ) is obtained, all other pertinent thermodynamic properties can 
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be readily evaluated by differentiating p(i uen ) with respect to the chemical potential in 
particular, the disorder-average mean particle density rS quen "> will be given by 

(9 

n (quen) _ p(quen) /g\ 

while the compressibility k T obeys 

,(quen) _ 1 ^ / R \ 

( n (quen)\ °V 

To close this section, we display the results corresponding to two "regular" cases: 
namely, when p = and p — 1, which will serve us in what follows as some benchmarks. 
In the p = all sites are decoupled, and one has the Langmuir results: 

pW = ^ln(l + z), n< ,on ) = -^-, and /T 1 A# m) = - (7) 
(3 v 7 1 + z 1 z K ' 

The "regular" case when p = 1 is a bit less trivial, but the solution can be still straight- 
forwardly obtained. In this case, we have 

/3 V 2 /' 1 + 4* - VTT4J' V ; 

and 

p-'k^ = - -. (9) 

T v / TT4i(i + 2z - v / TT4l) 

Note that in the p = 1 case (the completely catalytic system) the mean particle density 
tends to 1/2 as z — > oo (compared to n^ lan " ) — > 1 behavior observed for the Langmuir case), 
which means that the adsorbent undergoes "ordering" transition and particles distribution 
on the lattice becomes periodic revealing a spontaneous symmetry breaking between two 
sublattices. In the limit z — > oo the compressibility vanishes as k% e9) oc 1/ \fz compared 
to the Langmuir behavior k^ an ^ oc \ jz. 

3 Recursion relations for Zn((). 

Let us first introduce an auxiliary, constrained partition function of the form 

JV-2 

Z' N (() = Z N (()\ nN=1 = zJ2^= llni II ( X - & rn n i+1 ) (l - Civ-i n N -,), (10) 



{n;} i=l 



i.e. Z' N (() stands for the partition function of a system with fixed set ( = {Q} and fixed 
occupation of the site i — N, un — 1. Evidently, we have that 

z N (Q = z N . 1 (Q + z , N (Q. (ii) 

Next, considering two possible values of the occupation variable un-i, i-e. n^-i = and 
njv-i = 1, we find that Z' N (() can be expressed through Z N _ 2 (() and Z' N _ 1 (() as 

N-3 

Z' N (0 = z^>^=^ U(l-Qn t n t+l ) + 

{m} i=i 

N-3 

+ z 2 {l - Civ-i) ^^^ 2 ™ ! JJ (l - Ci n, ni+i) (l - Gv-2 njv-a) = 

{n<} i=l 

= zZ n ^(C) + z(1-Cn-i)Z' n ^(C) (12) 

Now, recursion in Eq.(ll) allows us to eliminate Z' N (() in Eq.(12). From Eq.(ll) we have 
Z'n(0 = Zn(() — Zn-i((), and consequently, we find from Eq.(12) that the unconstrained 
partition function Z N (() in Eq.(2) obeys the following recursion 

Z N {£) = (l + z(l - Ov-i))^v-i(C) + zCn-iZ N - 2 (C), (13) 

which is to be solved subject to evident initial conditions 

Z (C) = 1 and Z 1 (C) = l + z. (14) 

A conventional way (see, e.g. Ref. [20-22]) to study linear random three-term recursions 
is to reduce them to random maps by introducing the Ricatti variable of the form 

Rn(C) = (15) 

In terms of this variable Eq.(13) becomes 

R N (C) =(l + z(l - Cn-i)) + p Cn ~* , with R 1 (() = R 1 = l + z, (16) 

which represents a random homographic relation. Once Rn(C) is defined for arbitrary N, 
the partition function Z N (() can be readily determined as the product, 

N 

z N (C) = 11^(0, (17) 

1=1 
6 



and hence, the disorder- average logarithm of the partition function will be obtained as 



Before we proceed further on, we note that recursion schemes of quite a similar form 
have been discussed already in the literature in different contexts. In particular, two 
decades ago Derrida and Hilhorst [20] (see also Ref . [23] for a more general discussion) have 
shown that such recursions occur in the analysis of the Lyapunov exponent F(e) of the 
product of random 2x2 matrices of the form 



where Z{ are independent positive random numbers with a given probability distribution 
p(z). Equation (19) is related, for instance, to the disorder- average free energy of an Ising 
chain with nearest-neighbor interactions in a random magnetic field, and appears in the 
solution of a two-dimensional Ising model with row-wise random vertical interactions [24] , 
the role of e being played by the wavenumber 9. The recurence scheme in Eq.(16) emerges 
also in such an interesting context as the problem of enumeration of primitive words with 
random errors in the locally free and braid groups [25]. Some other examples of physical 
systems in which the recursion in Eq.(16) appears can be found in Ref. [22]. 

Further on, Derrida and Hilhorst [20] have demonstrated that F(e) can be expressed as 



N 




i=l 



(18) 




(19) 




(20) 



where R[ 



are defined through the recursion 



i?; = l + ^_i + ^_i(e 2 -l)/i2;_ 1 , with R[ 



= 1. 



(21) 



Moreover, they have shown that the model admits an exact solution when 



p(z) = (l-p)6(z)+p8(z-y), 



(22) 
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i.e. when similarly to the model under study, Zi are independent, random two-state vari- 
ables assuming only two values - y with probability p and with probability 1 — p. Suppos- 
ing that when % increases, a stationary probability distribution P(R') of the R\ independent 
of i exists [26], Derrida and Hilhorst [20] have found the following exact result: 

F(e) = pln(l + fo) -p(2-p)ln(l + 6^^-) + 

I -by 

°° / n — h \ N+l 

+ (l-P) 2 E^ ln (l + &(f^-) ), (23) 



by* 

N=l 



where 



1 + 



[i-y) 2 



1-1 + 4- 



e 2 y 



1/2- 



(24) 



2e*y L V {l-y) 2 J 

4 Disorder-average pressure. 

We turn now back to our recursion scheme in Eq.(16) and notice that setting 

R l {C) = {l + z)R' l , (25) 

and choosing 

-n ilan \ and e 2 = — ^- = ra< tan \ (26) 



1+z l+z 
makes the recursion schemes in Eqs.(16) and (21) identic! Consequently, the disorder- 
average pressure per site in our random catalytically-activated reaction/adsorption model 
can be expressed as 

p( q uen) = ]_ ^ + ^ + 1-F{e), (27) 

fj jj 

where F(e) is the Lyapunov exponent of the product of random 2x2 matrices in Eq.(19), 
in which e and Zi are defined by Eqs.(22) and (26). 

Note next that the first term on the right- hand- side of Eq.(27) is a trivial Langmuir 
result for the p = case (adsorption without reaction) which would entail n^" 6 ™-* — z/(l + 
z), Eq.(7). Hence, all non-trivial, disorder-induced behavior is embodied in the Lyapunov 
exponent F(e). 
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The disorder-averaged pressure per site for the random reaction/adsorption model un- 
der study can be thus readily obtained from Eqs.(23) and (24) by defining the parameters 
y and e as prescribed in Eq.(26). This yields the following exact result: 



where 



and 



pp(quen) = ln (^g _ (1 _ p ) _ ^ + 



^ N=l 

l + V / TT4^ 



<t>z 



VT+Tz-i 2 l 

uo = -. = z 6, = 1 — 

VT+Iz + 1 ' 2 2 



(28) 



(29) 



(30) 



Note that 2 obeys 4> z (4> z — 1) = z; hence, <p z= i = (\/5 + l)/2 is just the "golden mean". 



5 Asymptotic behavior of the disorder-average pres- 
sure, mean density and the compressibility. 

Consider first the asymptotic behavior of p(? uen ) in the small-z limit. To do this, it is 
expedient to use another representation of p(i uen ). After some straightforward calculations, 
one can cast 

p( q uen) in Eq.(28) into the form: 

^piemen) = (1 -P) ^p«jr n) 



V 



(31) 



n=0 



where T n denote natural logarithms of the Stieltjes-type continued fractions of the form 

/ \ 



T n = In 



1 + 



1 + 



1 + 



V 



1 + 



(32) 



1 + z / 



Note now that in the limit n — > oo, one has 



lim T n = ln(0 2 ) = In ( 



l + v / TT4i 



)• 



(33) 



i.e. T n is the n-th approximant of ln(0 2 ); hence, P^ uen ) can be thought of as the generating 
function of such approximants. Now, one finds that for z < 1 the sequence of approximants 
converges quickly to ln(</> 2 ); expanding the n-th approximant T n into the Taylor series in 
powers of z, one has that the first n terms of such an expansion coincide with the first n 
terms of the expansion of ln((f) z ), i.e. 



vrz; V 2 / 2v^F^ r(n+l) n V ; 

v n=l v ' 

Consequently, T n and T n -\ differ only by terms of order z n , which signifies that convergence 
is good. On the other hand, for z > 1 convergence becomes poor and one has to seek for 
a more suitable representation. As a matter of fact, already for z — 1 one has that in the 
limit n — ■> oo the approximant JF„ tends to hi(0i), i.e. the logarithm of the "golden mean", 
which is known as the irrational number worst approximated by rationals. 

In the small- z limit, we find then using an expansion in Eq.(34) that p(i uen ) follows 

ppiiuen) =z _ ^1 +p ^ z 2 + Q + 2p + p 2^3_ Q + 7 p + 4p2+p 3^4 + 0(z 5 ) (35) 

Consequently, in the small- 2; limit the mean density obeys: 



n (quen) = z 



(1 + 2p)z 2 + (l + Qp + 3p 2 ^ 3 - (l + 14p + 16p 2 + Ap^z 4 + C(^ 5 ), (36) 
while the compressibility k!^ uen ' 1 is given by 

Note that , the coefficients in the small- z expansion coincide with the coefficients in the 
expansions of p(' an ) and p( re f) when we set in Eq.(35) p — or p — 1. 

Now, we turn to the analysis of the large- ,2 behavior which is a bit more complex than 
z <C 1 case and requires understanding of the asymptotic behavior of the sum 

oo 

S = ^p Jv ln(l-(-l)V +2 ) (38) 

N=l 

entering Eq.(28). We note first that in this sum the behavior of the terms with odd and 
even N is quite different and we have to consider it separately. Let 

j oo 

p 



1 oo 

^ = -^p 27V ln(l + ^) (39) 



N=l 

10 



denote the contribution of the terms with odd N. Note that when z — > oo (i.e. uj — > 1) 
the sum S'odd tends to pm(2)/(l — p 2 ). The corrections to this limiting behavior can be 
defined as follows. Expanding In ^1 + uj 2N+1 j into the Taylor series in powers of uj and 
then, using the definition uj — 1 — l/(f> z and the binomial expansion, we construct a series 
in the inverse powers of (p z : 

s M = -L-mj)- 1 ^ 3 -^ 1 ■ 



1-p 2 w 2(l-p 2 ) 2 2 
l p(3 + 6p 2 -p 4 ) 1 1 p(15 + 10p 2 -p 4 ) 1 f4m 
8 (1-p 2 ) 3 2 24 (1 -p 2 ) 3 3 \<ft) 1 j 

Note that this expansion is only meaningful when <p z 3> (1 — p) -1 , (-2 3> (1 — p)~ 2 ), which 

signifies that p = 1 is a special point. 



Further on, plugging into the latter expansion the definition of <f> z , <p z = (1+y/l + 4z)/2, 
we obtain the following expansion in the inverse powers of the activity z: 

^ - r^ ln ^-2(i-p 2 ) 2 ^ + 

p (9-2p 2 + p 4 ) l p (3-4p 2 +p 4 ) 1 u v 

8 (1-p 2 ) 3 z 48 (1-p 2 ) 3 z 3 / 2 \z 2 ) 1 j 

Consider next the sum 

oo 

D_=£p^ln(l-c^+ 2 ), (42) 

N=l 

which represents the contribution of terms with even N . Note that in contrast to the 
behavior oiS dd, the sum in Eq.(42) diverges when z — > oo (uj — > 1). Since 1— w 2Ar+2 ~ 1— c<j 
when a; — > 1, we have that in this limit the leadin behavior of S even is described by 

2 

S even ~^j 2 \n{l-u). (43) 
To obtain several correction terms we make use of one of Gessel's expansions [27]: 

v ' k=i 

where gk(2N + 2) are the Dedekind-type sums of the form 

g k (2N + 2)= J2 7 — ~Tfc' ( 45 ) 

c 2JV + 2 =i,c^i (C - i j 
n 



where the summation extends over all ( being the (2N + 2)-th roots of unity (with ( — 1 
excluded). As shown in Ref.[27], the weights Qk{2N + 2) are polynomials in N of degree 
at most k with rational coefficients. Next, setting x = l/<p z in the expansion in Eq.(44), 
plugging it to Eq.(42) and performing summations over N, we find that S even can be written 
down as 

2 2 oo / lU 

S even = ln(^) + JL^ ln(2) + s P -J2 G *(p)^-, (46) 

where s p is an infinite series of the form 1 



^ = E^ ln ( Ar + 1 )' ( 48 ) 

N=l 

while Gk{p) are the generating functions of the polynomials gk(2N + 2): 

oo 

G k {p) = Y,9k{2N + 2)p 2N . (49) 



N=l 



Inserting next the definition of (f> z , we find the following explicit asymptotic expansion 



1 p 2 p 2 

Seven = o ln W + ~, o ln ( 2 ) + S P + 

2 1 — p z 1 — £) z 



p 2 (2-p 2 ) 1 p 2 (21-18p 2 + 5p 4 )l p 2 (2-p 2 ) 1 
(1 -p 2 ) 2 ^V2 + 24(1 -p 2 ) 3 I + 24(l-p 2 ) 2 2 3 / 2 



(?) < 50 > 



Finally combining the expansions in Eqs.(28), (41) and (50), we find the following large- ,2 
expansion for the disorder-averaged pressure p(i uen ) : 

v p(quen) = ih Hz) ~jr^ ln(2) + 

+ { —f^ s p +6(i + P y(i - p2) i + ° (^) > ( 51 ) 

which yields 



1+p 6 (1 + p) 2 (l — p 2 ) -2 \z 



1 Note that s p shows a non-analytic behavior when p — > 1. This function can be represented as 

ma-^-^E^^W)' (47) 



1 — p 2 I — p 2 n 

1 1 71=2 



where $(p 2 , n, 1) are the Lerch transcedents, $(p 2 , n, 1) = X^i=o(l + n P 21 - It is straightforward to find 
then that s p = —j^i ln(l — p 2 ) — + 0(\n(p)), where 7 is the Euler constant. 
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and 



^-l£,(<?«en) 



1 6 + 3p-p 3 1 i p(6 + 3p-p 3 ) i 
6(l+p)(l-p 2 )^ + 36(l+p) 2 (l-p 2 ) 2 ^ 




Note that asymptotic expansions in Eqs.(51), (52) and (53) are only meaningful for z ^> 
(1 — p)~ 2 and thus exclude the completely catalytic p — 1 case. At this special point p = 1 
we find from Eqs. (8) that the pressure and mean density exhibit a non-analytic dependence 
on 1/z: 



which differs substantially from the asymptotical behavior in the p < 1 case, Eqs. (51) 
and (52). This happens apparently because the bulk contribution to the disorder-average 
pressure in Eqs. (51) comes from the intervals devoid of the catalytic segments, in which 
reactions can not takes place and the mean density n ~ 1 in accordance with the Langmuir 
adsorption/desorption mechanism. Such intervals exist for any p strictly less than unity; 
their contribution vanishes only when p = 1. 

6 Conclusions. 

To conclude, in this paper we have presented an exact solution of a random catalytic 
reaction/adsorption model, appropriate to the situations with the catalytically-activated 
reactions on polymer chains containing randomly placed catalytist. More specifically, we 
have considered here the A + A — > reaction on a one-dimensional regular lattice which is 
brought in contact with a reservoir of A partilces. Some portion of the intersite intervals 
on the regular lattice was supposed to possess special " catalytic" properties such that they 
induce an immediate reaction A + A — > 0, as soon as two A particles land onto two vacant 
sites at the extremities of the catalytic segment, or an A particle lands onto a vacant 
site while the site at the other extremity of the catalytic segment is already occupied by 




(54) 



and 




(55) 
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another A particle. For quenched random distribution of the catalytic segments, we have 
determined exactly the disorder-averaged pressure per site and have shown that it can be 
represented as a sum of a Langmuir-type contribution and a reaction-induced term. The 
latter can be expressed as the Lyapunov exponent of a product of random two-by-two 
matrices, obtained by Derrida and Hilhorst [20]. Explicit asymptotic expansions for the 
mean particle density and the compressibility were also derived. 
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